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Abstract 


A fatigue crack closure model is developed that includes interactions 
between the three closure mechanisms most likely to occur at threshold; 
plasticity , roughness , and oxide. This model , herein referred to as the 
CROP model (for Closure , Roughness , Oxide , and Plasticity), also 
includes the effects of out-of-plane cracking and multi-axial loading. 

These features make the CROP closure model uniquely suited for ; but not 
limited to, threshold applications. Rough cracks are idealized here as 
two-dimensional sawtooths, whose geometry induces mixed-mode crack- 
tip stresses. Continuum mechanics and crack-tip dislocation concepts 
are combined to relate crack face displacements to crack-tip loads . 

Geometric criteria are used to determine closure loads from crack-face 
displacements. Finite element results , used to verify model predictions , 
provide critical information about the locations where crack closure 
occurs. 

Introduction 

Fatigue crack closure occurs when crack faces contact during cyclic loading. It is widely assumed that 
closed cracks are shielded from fatigue loads and no fatigue damage occurs while crack faces are in 
contact (ref. 1). Further, fatigue cracks are normally closed during a portion of the load cycle, so the 
applied AK is not a true crack-tip driving force. Instead fatigue crack growth (FCG) rate, dfl/dN, is 
related to an effective AK defined in Equations 1 and 2, where K c i is the stress intensity below which the 
crack is closed. Closure can even occur under tensile loading and has been used to explain load ratio 
effects on FCG (ref. 1). 


AKeff-Kmax K c) (1) 

du/dN = /(AK eff ) (2) 


Crack Closure Mechanisms 

Elber, the first to propose fatigue crack closure, suggested that crack face contact was due to residual 
plastic deformation in the wake of a propagating crack (ref. 1 ). This mechanism was coined plasticity- 
induced crack closure (PICC) as other closure mechanisms were suggested. Adams (ref. 2), 
Purshothaman and Tien (ref. 3), and Walker and Beevers (ref. 4) suggested contact of misaligned rough 
crack surfaces (roughness-induced crack closure, RICC) as a closure mechanism. The formation of a 
voluminous oxide layer along the crack faces (oxide-induced crack closure, OICC) was suggested as a 
closure mechanism by Walker and Beevers (ref. 4), Endo, et al. (ref. 5), and Suresh, et al. (ref. 6). For 
Paris regime FCG of most engineering alloys, PICC is the dominant closure mechanism ( i.e . RICC and 
OICC are second-order effects) and PICC models successfully predict closure loads and FCG rates (ref. 
7). Normally, RICC and OICC only become first-order effects as crack-face asperities and oxide layers 
become large relative to crack-tip displacements near the FCG threshold (ref. 8).* PICC is always a first- 
order effect and accounts for as much as 80% of the closure effects at threshold (ref. 9). Although 


* In cases of extreme crack path torturosity (ref. 10) and aggressive environment (ref. 5) RICC and OICC can be 
significant even in the Paris regime. 
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typically studied individually, it is likely that all three mechanisms (PICC, RICC, and OICC) are 
operative and interacting, especially near the FCG threshold. To accurately model closure at threshold, 
the combined influences of PICC, RICC, and OICC must be considered. 

Crack Closure Models 

Many models have been developed to predict fatigue crack closure behavior. However, most models 
are limited to a single closure mechanism and are unable to describe the interactions between multiple 
closure mechanisms likely to occur near the FCG threshold. Plasticity-based models adequately predict 
R-effects in the Paris regime where PICC dominates, but lack fidelity at threshold because RICC and 
OICC are neglected (ref. 9). RICC and OICC models have been developed to predict closure loads at 
threshold, but are limited because PICC is not considered, or they rely on empirical relations (which may 
be application specific) (refs. 1 1 and 12). Suresh and Ritchie were among the first to model rough cracks 
as a two-dimensional sawtooth wave (ref. 1 1). This rough crack geometry produces mixed-mode crack- 
tip stresses, and both opening and sliding mode crack-face displacements. As modeled by Suresh and 
Ritchie, mode II displacements do not return to zero upon unloading, creating misalignment between 
mating crack surfaces. Wang, et al. (ref. 12) developed a similar model, where mode II displacements 
were calculated in terms of dislocations emitted from a single slip plane coinciding with the crack-tip 
direction. Because the Wang model only allows slip to occur along a single plane, only sliding mode 
displacements occur, although mixed-mode crack-tip stresses are expected. However, neither model 
considers crack-wake plasticity or oxides. 

Llorca (ref. 13), Garcia and Seheitoglu (refs. 14 and 15), Chen and Lawrence (ref. 16), and 
Ravichandran (refs. 17 and 18) have proposed more complicated RICC models (in terms of computation 
requirements). Llorca developed a finite difference model, which also idealized the crack profile as a 
sawtooth wave. The most important variables influencing closure levels for this model were crack tilt 
angle and average direction of crack growth (considering out-of-plane cracking). Furthermore, RICC was 
discontinuous; for instance, closure was greatest as the crack was changing direction. This observation is 
significant because it suggests closure is more likely to occur at specific crack-wake features (i.e., at 
sawtooth peaks and valleys). Garcia and Seheitoglu modeled rough crack asperities as three-dimensional 
features randomly distributed along the crack surface. Statistical means were used to evaluate asperity 
heights and radii. The crack surface was divided into multiple strips, and each strip was treated as a 
rough surface contact problem. Chen and Lawrence investigated the combined effects of RICC and PICC 
by altering a strip yield model, developed by Newman (ref. 19), for a sawtooth crack. The Chen and 
Lawrence model is among the first to describe the interaction of closure mechanisms (RICC and PICC). 
The Ravichandran model also used a sawtooth profile and predicted closure levels as a function of 
microstructural scale (i.e., grain size), dislocation density, and yield stress. As with the Wang model, 
sliding mode displacements were a result of crack-tip dislocation emission in the Ravichandran model. 

Newly developed numerical models (refs. 13-18) are more complex and sophisticated than the 
analytical models of Suresh and Ritchie (ref. 11) and Wang, et al. (ref. 12); some are able to consider 
multiple closure mechanisms. Despite the added complexity, no model adequately predicts the combined 
effects of PICC, RICC, and OICC, both near the FCG threshold and in the Paris regime. Ideally, a 
closure model would be complex enough to consider the combined effects of PICC, RICC, and OICC (the 
closure mechanisms considered most likely at threshold), yet be computationally simple and analytical in 
nature. 
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Objectives 

The objective of this paper is to develop an analytical fatigue crack closure model that predicts and 
describes contributions and interactions of PICC, RICC, and OICC; considered the most likely closure 
mechanisms at threshold. The closure model developed in this paper will be called the CROP model for 
Closure, Roughness, Oxide, and Plasticity. Special considerations will be made for near-threshold FCG 
scenarios, but this model will not be limited to low AK. The crack geometry is idealized as a sawtooth 
wave, but modified to allow out-of-plane cracking. To develop this model, it is necessary to characterize 
crack-face displacements of rough cracks, and include crack-wake plasticity and oxide. Cracks might 
kink on a macroscale due to mixed-mode loading or geometry (ref. 20). However, it is the development 
of kinks on a microscale (/. e. roughness), due to microstructural features (ref. 21), that are of primary 
interest in this paper. Most existing closure models include a single closure mechanism and work well for 
either Paris or threshold FCG, but not both. By including multiple closure mechanisms, it is hoped that 
the CROP model will work for a broad range of FCG scenarios. Finally, if supported by experimental 
data (see ref. 22), the CROP model can be used to study interactions between PICC, RICC, and OICC. 

Model Development 

Step-by-step development of the CROP closure model is presented in this section. In addition to 
modeling PICC, RICC, and OICC, out-of-plane cracking and a limited amount of (applied) mixed-mode 
loading are considered. f Inclusion of each mechanism is discussed separately. Crack-face displacements 
are characterized and compared with crack-wake effects (i.e., plasticity, roughness, and oxide) to compute 
closure loads. 

Consideration of RICC and out-of-plane cracking 

Rough crack surfaces are modeled here as a sawtooth wave (i.e., a series of two-dimensional, rigid, 
triangular asperities, of equal size and shape) similar to the geometry of existing closure models (refs. 1 1- 
13, 16, 17). A schematic of this sawtooth crack profile, modified to consider out-of-plane cracking (a 
macroscopic deviation from global coordinates by an angle, y), is shown in Figure 1.* Characteristic 
roughness parameters (asperity length, g, and asperity angle, a) should be obtained from real crack 
surfaces, or related to microstructural features, such as grain size. The sawtooth geometry is a reasonable 
approximation for most rough crack surfaces, and simplifies the geometric conditions of crack face 
contact. Three coordinate systems are required to describe the crack geometry in Figure 1 . First, a global 
coordinate frame is used to describe loading or specimen geometry, and from which out-of-plane crack 
deflection is measured. The global coordinate frame is labeled by X and Y in the figure. Next, a local 
coordinate frame is needed. The x-axis of the local coordinate frame coincides with the nominal crack 
direction (average direction of FCG from a macroscopic point of view). This coordinate frame is labeled 
in Figure 1 by x and y. Finally, a crack-tip coordinate frame is needed, whose x-axis coincides with the 
orientation of the crack-tip on a microscale. As shown in Figure 1, the crack-tip coordinate frame is 
labeled by x' and y'. When considering closure of rough cracks, it is convenient to describe crack face 
displacements using both local and crack-tip coordinates. 


f The CROP model is limited to cases where mode I loading is dominant. Problems arise when sliding-mode crack 
closure (SMCC) occurs, a result of mode II loading, as described in detail elsewhere (ref. 23). 

* ASTM standards permit neglecting out-of-plane cracking at angles less than 10°, and FCG data for angles up to 20° 
is considered valid, but out-of-plane cracking must be reported (ref. 24). 
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Figure 1. A schematic of the idealized crack configuration used by the CROP model; a two-dimensional sawtooth 

subject to out-of-plane cracking. 


Stress intensity factors 

Crack-tip stress-intensity factors for out-of-plane rough cracks are approximated using analytical 
stress solutions for two crack configurations. Solutions of angled straight cracks and angled, infinitesimal 
kinks at the tip of an otherwise straight cracks are used to predict the effects of out-of-plane cracking and 
crack roughness, respectively. The methodology used to determine crack-tip stresses for rough out-of- 
plane cracks is shown schematically in Figure 2. First, the analytical solution for an angled straight crack 
is used to predict the effects of out-of-plane cracking. This crack configuration is shown in Figure 2a 
under general planar loading, i.e., with non- zero cty, <Jx, and txy. Crack-tip stress-intensity factors are 
obtained for this configuration by transforming the applied loads (global coordinates) into local 
coordinates. These transformed stresses are then used to calculate the stress-intensity factors shown in 
Equations 3 and 4, where Ki g = c Y Vrca, Kn g = XxyVrta, and K 0g = Gx'lna. The stress intensity factor Ko g 
accounts for non-singular stresses parallel to cracks that contribute to the singularity upon meandering 
(refs. 25 and 26). Note that 2 a is defined as the crack length projected onto the plane perpendicular to o Y ; 
and is used because this is the dimension normally measured during fatigue crack growth (FCG) tests. 
The actual crack length is 2a/cos (y). 


^ i Local = K lg Vcos 3 (y) -2K Ug sin(yX/cos(y) + K 


Og 


sin 2 (y) 


\/ c °s(y) 


(3) 


^ u Local =K Ig sin(yycos(Y)-K 


ng 


cos 2 (y) 
Vcos(y) 


-KogSin(y)Vcos(y) 


(4) 


The stress intensity factors for a straight crack with an infinitesimal kink are used to approximate 
stress-intensity factors for rough (sawtooth) cracks (see Figure 2b) (ref. 27). These relations are given as 
Equations 5 and 6, but are valid only if the kink length, 5 a, is small compared to the K-dominated region 
at the crack-tip (ref. 28). It will be shown later that Equations 5 and 6 are also valid for sawtooth cracks. 
Note that local stress intensity factors of Equations 3 and 4 are used in Equations 5 and 6 allowing 
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(a) Angled straight crack 


Kn Local 



Figure 2. Schematics to describe analytical stress-intensity factor solutions for two crack configurations - an angled 

straight crack and a straight crack with a kinked crack tip. 


simultaneous consideration of crack roughness and out-of-plane cracking. As schematically shown in 
Figure 2, the angled straight crack solutions are used first to determine the macroscopic effects of out-of- 
plane cracking (Equations 3 and 4), then, the microscopic crack roughness effects are evaluated 
(Equations 5 and 6). The stress intensity factors are expressed in terms of globally applied stress-intensity 
factors because these quantities are normally used to describe FCG tests. 
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Elastic crack face displacements 

Elastic crack-face displacements are functions of the crack-tip stress-intensity factors determined in 
the previous section. Because all rough cracks are subjected to a mixed-mode stress state, both opening 
and sliding mode crack-face displacements (mode I and mode II, respectively) are expected. Sliding 
mode displacements create misalignment of crack-face asperities, i.e., RICC. Elastic crack face 
displacements for straight cracks under plane-strain conditions and mixed-mode planar loading are shown 
in Equation 7, where E is the elastic modulus, v is Poisson's ratio, and r is the distance behind the crack 
tip (see Figure 1) (ref. 26). Note that / = r / cos(a) and the total relative displacements between the two 
crack faces are 2u and 2v'. Primed quantities indicate use of crack-tip coordinates and r is defined with 
respect to the local coordinate frame. Plane-strain conditions were adopted because the CROP model is 
intended for threshold applications; for typical FCG specimens near threshold, most of the crack front is 
under plane-strain conditions. It will be shown later that Equation 7 is valid in the crack-tip vicinity for 
sawtooth cracks. The crack-face displacements of Equation 7 are expressed in crack-tip coordinates, but 
can be transformed to local coordinates using Equation 8. Local coordinates will be convenient when 
considering crack closure. 
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Consideration of crack-tip plasticity and PICC 

As a cracked body is loaded, a plastic zone forms at the crack-tip (refs. 29 and 30). Crack-tip 
plasticity for mixed-mode loading will produce a crack-tip opening displacement (CTOD) and a crack-tip 
sliding displacement (CTSD). § Existing crack-tip slip models for mode I straight cracks (refs. 3 1 and 32) 
consider all crack-tip deformations as the result of dislocation motion along planes that coincide with the 
maximum singular shear stress. Although these models were developed for mode I loading, this approach 
is generalized by determining the orientation of extreme values of the singular crack-tip shear stresses for 
mixed-mode loading are modified to estimate plastic displacements of a mixed-mode crack tip. On a 
microscale, plasticity is the result of dislocation glide originating near the crack tip (ref. 33). Near the 
crack-tip, dislocation pairs of opposite sign are created. Dislocations of one sign are attracted to the crack 
tip, where they are annihilated, producing plastic displacements. Dislocations of the opposite sign are 
driven away from the crack-tip. Ultimately, these dislocations will either (1) stop at a barrier, such as a 
particle interface or grain boundary (refs. 34-36), or (2) go far from the crack-tip where a reduction in 
dislocation driving force (stress) occurs (refs. 33 and 37). For the model developed herein, slip is 
assumed to occur by edge dislocations that glide freely along slip paths until reaching the edge of the 
plastic zone, which acts as a Peierls-Nabarro energy barrier; here stresses become too small for further 
dislocation movement. As in previous crack-tip dislocation emission models, dislocations are assumed to 
travel along planes originating at the crack-tip.** The plane-strain plastic zone radius is defined in 
Equation 9, as a function of angle from the crack plane, 0 (ref. 39). Equation 9 is based on linear-elastic 
material behavior and assumes no stress redistribution due to crack-tip plasticity; this radius is that of the 
elastic solution whose vonMises equivalent stress corresponds to yielding. The increase in plastic zone 
due to stress redistribution was accounted for by increasing the plastic zone of Equation 9 by a factor of 
1.828 (see Equation 10). This constant was chosen such that the CTOD (which will be calculated using 
Equation 9) for a straight crack matched the theoretical value (K 2 / E a 0 ). This relation is used to 
determine the monotonic plastic zone (MPZ) associated with K max . n 

It is recognized that slip occurs along preferred slip planes that depend on microstructure and local 
shear stress (acting on these planes). However, it is assumed that slip planes (on average) coincide with 
planes of extrema in elastic crack-tip shear stress, This combination of dislocation emission and 
continuum mechanics is appropriate when crack-tip deformation occurs in multiple grains along the crack 
front. The elastic, singular crack-tip stresses for mixed-mode loading under plane-strain conditions are 
given in Equation 1 1 (ref. 26). The directions of shear extrema, 0,, satisfy Equation 12 (which is obtained 


§ In general, CTOD and CTSD may include both elastic and plastic components. However, Equation 7 indicates no 
elastic displacement exists at the crack tip. For this model CTOD and CTSD describe only plastic crack-tip 
displacements. 

** Slip is described as occurring on planes whose projection into the two-dimensional model is a line. However, the 
phrase "slip plane" is used to be consistent with existing literature (ref. 38). 
n The MPZ is the crack-tip region where plasticity occurs during loading, and is associated with K nux . The cyclic 
plastic zone (CPZ), a subset of the MPZ where reverse plasticity occurs during unloading, is associated with AK. 
The CPZ will be discussed later. 
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by setting the derivative of x^, with respect to 0 to zero). This equation is transcendental and must be 
solved numerically, i.e. no closed form solution exists for 9j. Equation 12 has three solutions except for 
straight cracks under pure mode I loading, which have only two solutions. A plot of the singular shear 
stress coefficient (x^ multiplied by \2jir) versus angle from the crack-tip direction is shown in Figure 3a 
for a crack configuration described by a = 45°, y = 10°, K Ig = 1, and K IIg = 0. The three local extrema are 
indicated by vertical dotted lines and labeled in Figure 3a. A schematic of the corresponding plastic zone 
(calculated by Equation 10) is shown in Figure 3b. The slip planes indicated in Figure 3a are shown in 
Figure 3b as solid lines emanating from the crack tip. Although the extrema of crack-tip shear do not 
exactly coincide with extrema in plastic zone radius, the two differ by less than 30°. 



The dislocation distribution along individual slip planes has been solved assuming a constant shear 
stress in the plastic zone (e.g., assuming elastic-perfectly plastic behavior at the crack-tip, use x = 0.5 <J 0 ) 
and no interaction between slip planes. The total slip along a slip plane of length r p is given by Equation 
13, where E is the elastic modulus, v is Poisson's ratio, n is the total number of edge dislocations, and b is 
the magnitude of the Burger's vector (ref. 40). 


B = nb = 



(13) 


Using this relation, plastic crack-tip deformations (CTOD' and CTSD') are computed by summing the 
vector components of each slip plane (in crack-tip coordinates) as shown in Equations 14 and 15. (Again, 
primed quantities correspond to crack-tip coordinates.) The symbol I indicates summation over the range 
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0 (degrees) 


(a) Singular crack-tip shear stress variation 



(b) Shape of corresponding plastic zone 



Figure 3. The effect of mixed-mode loads on crack-tip stresses and the crack-tip plastic zone. 


of the index ‘i\ which typically ranges from one to three (the total number of slip planes). ** Conversion 
to local coordinates involves the same transformation given in Equation 8. 


CTOD'= J]B i sin(e i ) 
CTSD' = ^ BjCOS^j ) 


04) 

(15) 


The cyclic plastic zone (CPZ) is a subset of the MPZ, where reverse plasticity occurs during 
unloading. The size of the CPZ is proportional to AK 2 , and the shape can be calculated using Equations 9 
and 10, by inserting one-half of the cyclic stress intensity factors in place of the maximum values (e.g., Ki 
Local is replaced by Vi AKi Local)- Reversed plasticity in the CPZ leads to a reduction in residual CTOD and 
CTSD. For modeling purposes, reversed slip is assumed to occur along the same slip planes as for 
loading, but using the CPZ boundary as maximum slip length and replacing x in Equation 13 with -x. 
Because the MPZ and CPZ radii are proportional to (K max ) 2 and VifAK) 2 , respectively, CTOD and CTSD 
are also a function of K max and AK, as in Equation 16. The subscript ‘max’ corresponds to maximum 
values of that quantity, which occur at maximum load. §§ 

JCTODl 
[CTSDJ 


1 - 


AK 


\ 2 


V^K-max J 


|CTOD„ 

[CTSD„ 




Although RICC and OICC are most significant at FCG threshold, PICC is still a major contributor and 
must be considered (refs. 9 and 41). Thus, a complete discussion of crack plasticity must include crack- 
wake deformation. A PICC model by Budiansky and Hutchinson (ref. 42) relates crack-wake plasticity to 


** Only two slip planes exist for a straight crack (i.e., a = 0), so T ranges from 1 to 2 in this case. 

55 This is only true in cases of proportional loading - that is, the ratio Kn/Ki remains constant during the load cycle - 
which is an assumption of the CROP model. 
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the CTOD. A value less than CTOD™* (the CTOD at maximum load) is used to account for reversed 
crack-tip plasticity during unloading. For R = 0, Budiansky and Hutchinson calculate the residual plastic 
wake as approximately 86% of CTOD max , or £ = 0.86. The normalized crack- wake plasticity, £ is defined 
in Equation 17, where CTOD wakc is the total residual plastic wake (half associated with each crack face). 
The relation between CTOD max and CTOD wake is shown schematically in Figure 4a. Here, a straight crack 
is shown under load with a blunted crack tip. Residual crack-wake plasticity is indicated as shaded 
regions along each crack face. Although plastic deformation partially fills the crack mouth in this 
schematic, the crack is shown open. Using results from Budiansky and Hutchinson (open circular 
symbols), a plot of £ versus applied R is shown in Figure 4b along with a least squares fitted (4 th order) 
polynomial. This fitted polynomial is expressed symbolically in Equation 18, and is used to estimate ^ as 
a function of R. Using Equation 18 instead of the entire procedure developed by Budiansky and 
Hutchinson simplifies evaluation of the CROP closure model. 


CTOD wake 

CTOD_ 


(17) 


£ = 0.856 1 + 0.0205 • R+ 0.1438* R 2 + 0.2802 R 3 - 0.3007 • R 4 


(18) 



(a) Crack tip at maximum load (b) Variation of £ with R (ref. 42) 

Figure 4. The relation between crack-wake and crack-tip plasticity. 


Consideration of OICC 


A crack mouth oxide layer is modeled as an incompressible wedge of uniform thickness (t) along both 
crack surfaces, similar to a model developed by Suresh (ref. 6). The voluminous oxide layer partially fills 
the crack mouth and increases closure levels. An oxide layer of finite thickness will produce higher 
closure levels than at higher values of K max because the available crack opening is small at low values of 
K max . Consequently, OICC is thought to be most significant near the FCG threshold (ref. 43). Naturally 
forming oxide layers on fresh aluminum surfaces have been measured as 10 A thick in dry air (ref. 44), 
and between 50 A and 150 A in humid air (ref. 45), although considerable variation with environment and 
temperature was noted. Fretting contact of rough crack surfaces (/.<?., RICC) can create a large amount of 
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oxide debris, a result of repeated removal and regeneration of crack-mouth oxide (ref. 43). Large 
amounts of oxide debris may also form in aggressive (oxide producing) environments (ref. 5), where the 
effective oxide thickness can be much larger: up to 1000 A in moist air for aluminum alloys (ref. 43) and 
3-5 pm for steels (ref. 46). Although this model assumes a uniform oxide layer on the crack faces, 
modification can be made to simulate oxide debris (non-uniform). For example, oxide debris generated 
by fretting asperities may accumulate behind the crack tip. Oxide debris normally is thicker than an 
undisturbed oxide layer, expected at the crack tip. In this case, closure can be determined separately, at 
the location of the oxide debris and at the crack tip, using different values of oxide thickness. The crack 
face location with the largest closure load will contact first and dominate crack-tip behavior. 

Numerical Verification of Analytical Relations 

Finite-element analyses were performed to validate crack-tip stress-intensity factor and displacement 
solutions (Equations 3-8) for sawtooth cracks (Figure 1) using FRANC2D (ref. 47). The asperity length, 
g , and asperity angle, a, were varied for the sawtooth configuration. Only mode I loading was considered 
(Kng = 0) with no crack meandering (y = 0) because this is the most common configuration for laboratory 
testing. Stress-intensity factors were extracted from the elastic solutions using the J-integral technique 
(ref. 48). A typical deformed (grossly exaggerated) mesh used for these analyses is shown in Figure 5 for 
a = 45° and g/a = 0.05. Taking advantage of symmetry, only half of the specimen is analyzed (using 
appropriate boundary conditions); the entire finite-element mesh is shown in Figure 5a. A rosette of 
singular quarter-point elements is used around the crack tip. The rest of the specimen is meshed with 
three- or four-sided quadratic elements. As seen in Figure 5b, the mesh is refined near the crack tip. At 
least four elements were used between the crack tip and the asperity nearest the crack tip, regardless of 
the asperity length, g. Further refinement of the crack-tip region changed the results by less than 2%. 
The sawtooth wave was extended to form at least 10 complete asperities behind the crack tip, or to the 



(a) Entire mesh (b) Refined mesh near crack tip 

Figure 5. A typical (deformed) finite-element mesh is shown for a = 45° and g/a = 0.05. 
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specimen centerline. For convenience, the remaining crack profile was modeled as a straight line; further 
extension of asperities along the crack surface changed results near the crack tip (0 < r/g < 1 ) by less than 
1 %. 

Stress-intensity factors obtained from finite-element analyses are plotted in Figure 6 for asperity 
angles (a ranging from 0° to 45°, and for normalized asperity length (g/a) ranging from 0.2 to 0.002. 
Mode I (K]) results are shown in Figure 6a, and mode II (Kn) results are plotted in Figure 6b. Finite- 
element results are shown as open symbols, while analytical solutions of Equations 5 and 6 ( g/a = 0) are 
shown as closed symbols. Refer to the legend in the upper left corner of the figure for symbols 
corresponding to specific values of a. This crack configuration is shown schematically in the upper right 
comer of the figure. As the normalized asperity length (g/a) approaches zero, the finite-element sawtooth 
crack results converge to the analytical kinked-crack solutions. This result validates the use of Equations 
5 and 6 for g/a < 0.05. Similar agreement between finite-element and analytical solutions were obtained 
for mode II loading, but results are not presented here (ref. 23). 



(a) Mode I results (b) Mode II results 

Figure 6. Stress-intensity factors for sawtooth cracks (under mode I applied loading) for a = 0°, 15°, 30°, and 45°. 


Numerical elastic crack-face displacements are plotted against r/g (normalized distance behind the 
crack tip) in Figure 7 for asperity angles of 30° and 45°. Note that all quantities are presented in 
dimensionless form where appropriate. Dimensionless values for crack-face displacements (u and v*) are 
defined as follows. 
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u 


(19) 


u E 
vE 

kJi-^VF 


( 20 ) 


Finite-element results for two values of gla (0.1 and 0.01) are shown as open symbols (circular and 
square, respectively). Dashed and solid lines correspond to analytical elastic solutions of sliding and 
opening mode displacements predicted by Equation 7 (transformed to local coordinates), respectively. 
The dotted line is the (mode I) displacement solution for a straight crack (i.e., a = 0) under the same 
mode I loading. As seen in Figure 7, mode II displacements (u) match the analytical solution well, but 
only between the crack tip and the first asperity behind the crack tip (0 < rig < 1). Further behind the 
crack tip ( rig > 1), mode II displacements exhibit an oscillatory decay about zero with extreme values 
corresponding to changes in crack direction, i.e., peaks and valleys of the sawtooth. The results for a = 
30 ° and a = 45 0 are similar in character, but the a = 45° results show a greater mode II displacement am 
plitude. The mode I displacements (v‘) also match the analytical solutions, but far from the crack tip the 
numerical solutions approach those expected of a straight crack. The finite element results in Figure 7 
suggest mode II elastic displacements are greatest at the asperity nearest the crack tip, but mode I 



rig rig 

(a) a = 30° (b) a = 45° 


Figure 7. Elastic crack-face displacements of sawtooth cracks (under mode I loading) from finite-element analyses 

and analytical model solutions. 
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displacements are smaller here than at any other asperity in the crack wake. Consequently, RICC is most 
likely to occur at the asperity nearest the crack tip. Where crack roughness effects are important only this 
first asperity needs to be considered; asperities further behind the crack tip can be neglected. This 
observation is consistent with finite difference results of Llorca, who reported mode II displacements 
were greater at asperity peaks and valleys than intermediate points (ref. 13). However, Llorca did not 
indicate that the asperity nearest the crack tip was the location where RICC was most likely. 

Evaluation of Closure Levels 

For rough cracks under applied mode I loading, finite element results reveal that crack-face contact 
might first occur at the asperity nearest the crack-tip. However, initial contact is expected at the crack tip 
for flat cracks (no roughness) because the elastic displacements are zero (ref. 49). In general, crack 
closure can first occur either at the crack tip (called tip contact) or at the asperity nearest the crack tip 
(called asperity contact). Simple geometric relationships are established here to determine crack closure 
levels for both types of crack face contact. A schematic of a sawtooth crack with an oxide layer is shown 
in Figure 8, where small boxes enclose the crack-tip and first asperity (a distance g behind the crack tip). 
Open arrows indicate high magnification schematics of these locations at the bottom of the figure. For 
asperity contact (left side of figure) the total opening and sliding displacements are labeled as (2v + 
CTOD) and (2 u + CTSD), respectively, and oxide layer thickness is indicated ( t ) as the shaded region. 
For convenience, local coordinates are used to describe asperity contact. For tip contact (right side of 
figure), only plastic displacements occur at the crack tip, so only CTOD' and t are labeled in the lower 
right corner of Figure 8. Crack-tip coordinates are used to describe tip contact. The geometric relations 
satisfied for asperity contact and tip contact are shown in Equations 21 and 22, respectively. The £ terms 
in these equations account for residual crack wake plasticity. In Equation 21, the absolute value of the 
mode II quantity (2w c i + CTSD) is used because sliding mode displacements in either direction promotes 
closure. Quantities with the subscript 'cf correspond to values evaluated at crack face contact. Because 
elastic displacements are proportional to K and plastic displacements depend on (K max ) 2 and (AK) 2 , all 



(a) Asperity contact (b) Tip contact 

Figure 8. A schematic of a rough sawtooth crack (top of figure) to highlight the two locations where closure may 
first occur; at the asperity nearest the crack tip or at the crack tip. 
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displacements can be expressed in terms of load ratio at closure (R d ), as in Equations 23 and 24. 
Combining Equations 21-24 allows the closure level, R d = K d / K max , to be determined at both critical 
locations. ” " Contact at either location will affect the crack-tip driving force. Therefore, the greater of the 
two closure levels calculated from Equations 21 and 22 is used as the R d affecting crack-tip behavior. 


(2v d + CTOD tl )cos(a) - |(2w d 


+ CTSD c , )sin(a)j = 2t + q— — — j- 


( 21 ) 


CTOD tl =2r + i;(CTOD m j' 


( 22 ) 



|CTOD d 

[CTSD d 



|CTOD n 

CTSD„ 


(23) 


(24) 


Summary 

A fatigue crack closure model (called the CROP model) has been developed that includes plasticity-, 
roughness-, and oxide-induced crack closure mechanisms (PICC, RICC, and OICC, respectively). The 
CROP model is uniquely suited for a wide range of AK, including threshold, because it accounts for the 
three closure mechanisms most likely at FCG threshold. Finite-element results indicate that RICC is most 
likely to occur first at the asperity nearest the crack tip because mode II displacements are greater and 
mode I displacements are less than at any other crack-wake asperity. However, PICC models indicate 
closure first occurs at the crack tip where crack roughness is not a factor. The CROP model considers 
both the crack tip, and the asperity nearest the crack tip, as potential locations where closure first occurs 
(called tip contact and asperity contact, respectively). Non-linear interactions between PICC, RICC, and 
OICC predicted by the CROP model are verified with experimental data in reference 22. 
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